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Can the 33 s pulsations observed from AE Aquarii be explained 
in terms of accretion onto the white dwarf surface ? 
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Abstract. The 33 s pulsing component observed from AE Aqr is frequently assumed to be the result of accretion 
onto the surface of a rotating white dwarf. The validity of this assumption is discussed in the present paper. I show 
that under the conditions of interest the white dwarf is in the state of supersonic propeller and the efficiency of 
plasma penetration into its magnetosphere is ^ 0.1%. This is too small to explain the observed luminosity of the 
pulsing component. Moreover, I find that for the currently established value of the angle between the magnetic 
and rotational axes of the white dwarf, the material entering the magnetosphere at the boundary can never reach 
its surface by fiowing along the magnetic field lines. I conclude that the assumption about the direct accretion 
onto the surface of the white dwarf in AE Aqr contradicts the observational data obtained in the optical/UV and 
thus cannot be accepted. 
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1. Introduction 

AE Aquarii is a non-eclipsing close binary system with 
an orbital period Porb ~ 9.88 hr and orbital eccentricity 
e « 0.02. It is situated at the distance of ^ (100 ± 30) pc. 
The normal component (secondary) is a K3-K5 main se- 
quence red dwarf. The primary is a magnetized white 
dwarf rotating with a period of 33 s. The inclination angle 

o o 

of the system is limited to 50 < i < 70 and the mass ra- 
tio isq = M 2/M1 = 0.77 ± 0.03 (for references see Table 1 
in Ikhsanov pOOOl ). 

AE Aqr emits detectable radiation in almost all parts 
of the spectrum. It was observed as a powerful non- 
thermal flaring source in the radio (Bastian et al 
and gamma-rays (Bowden et al. |1992| : 
1994). The optical, UV and X-ray radiation of the sys 



1988D 
Meintjes et al. 



(T ~ 26000 K) spots with a projected area of~ 4 10^^ cm^, 
which are situated on the white dwarf surface. Placing 
these spots at the magnetic pole regions of the primary, 
Eracleous et al. ( 1994 ) have estimated the value of the an- 
gle between the rotational and magnetic axes of the white 

o o 

dwarf to be 75 ^ (3 <, 77 . The remaining light comes 
from a highly variable extended source which manifests 
itself in the blue/UV continuum, the optical/UV broad 
single-peaked emission lines and the non-pulsing X-ray 
component. This source is spread around the magneto- 
sphere of the primary and is associated with the mass- 
exchange process in the system (for a detailed syst em de - 
scription see Eracleous & Horne 1996 and Ikhsanov 2000). 



tem is predominantly thermal and is coming from at 
least three different sources. The visual light is dominated 
up t o 95%) by the secondary (Bruch 1991 Welsh et al. 
1995 ) . The contribution of the primary is observed mainly 



in the form of 33 s (and 16.5 s) coherent oscillatio ns de- 
tectable in the optica l, UV and X-rays (Patterson 1979 ; 
Patterson et al. 19801) . Study of the optical- UV spectrum 
and profiles of the pulsations (Eracleous et al. 1994) has 
shown that the pulsing radiation is coming from two hot 



On the grounds of the properties observed in the 
optical, AEAqr was classified as a nova-like star be- 
longing to the DQ Herculis subclass of magnetic cata- 
clysmic variables. Following this notion, it has been sug- 
gested that the primary is a magnetized accretion-powered 
white dwarf undergoing disk accretion (Patterson 1979). 
Correspondingly, the flaring activity of the system was in- 
terpreted in terms of a variable mass accretion rate onto 
the primary surface (van Paradijs et al 



1989) 
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However extensive investigations of the system over the 
last five years have clearly shown that AE Aqr does not fit 
this model. First, studies of th e 33 s pulsations in the opti- 
cal/U V (Eracleous et al. 1994) and X-rays (Reinsch et al. 

1999| ) have 



1995; Clayton & Osborne 1995; Choi et al 
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shown no correlation between their ampUtudes and the 
flaring of the system. This allowed Eracleous et al. (1994) 
to conclude that flares are not related to the process of de- 
positing material onto the primary surface. Second, anal- 
ysis of the Ha Doppler tomogram of AEAqr has shown 
no evidence of a developed Keplerian accretion disk in the 
system. Instead, it has shown that the material inflowing 
through the LI point into the Roche lobe of the primary 
is then streaming out from the system with an average ve- 
locity - SOOkms-i ( Wynn e t al. |1997t Welsh et al. |1998D . 
Finally, de Jager et al. (1994) have reported the mean spin- 



down rate of the white dwarf Pg = 5.64 x 10^ s s~ which 
implies the spindown power of Lgd — I^Cl ^ lO'^^ergs^^. 
Lsd exceeds the observed UV and X-ray luminosities of 
the system by a factor of 120 and its bolometric luminos- 
ity by a factor of more than fiveQ. This indicates that the 
spindown power dominates the energy budget of the sys- 
tem and raises the question of the nature of the spindown 
torque, which would appear to be much larger than any 
inferred accretion torque. 

Four alternative answers to this question and, corre- 
spondingly, four alternative theoretical models of the sys- 
tem are currently under dis cussio n: (i) the magnetic pro- 
peller model (Wynn et al. 1997 ), (ii) th e spin-powered 
white dwarf pulsar model (Ikhsanov 199S ), (iii) the gravi- 
tational wave emitter model (Choi & Yi 2000 ) and (iv) the 
differentially rotating white dwarf model (Geroyannis 
2001). Though these models essentially differ from each 
other in many aspects they have an important common 
point. Namely, the mass-exchange picture in AEAqr is 
interpreted in terms of the propeller action by the white 
dwarf. This implies that the magnetic field of the white 
dwarf is strong enough for the magnetospheric radius. 



R 



2/7 



to be larger than its corotational radius. 



Rc 



1-410^ Mo.i'i^'f cm. 



(1) 



(2) 



Here ^, Mq.s = Mi/O.SM© and P33 = Ps/33s are the 
magnetic dipole moment, the mass and the spin period 
of the white dwarf, respectively. M is the mass accretion 
rate onto the primary magnetosphere and 0.5 ^ ^ 1 
is the parameter which accounts for the geometry of the 
accretion flow. Under these conditions the white dwarf is 
in the centrifugal inhibition regime and is able to eject the 
material out of the system due to the propeller action. 

One of the important questions, which remains to be 
explained within this approach, is the origin of the hot 
spots on the white dwarf surface. According to the canon- 
ical propeller model (e.g. Davies et al. 197S| ) no accretion 
occurs onto the surface of a star in the state of propeller. 
This indicates that the direct accretion onto the surface 
of the white dwarf in AEAqr cannot be responsible for 
the origin of the hot polar caps and hence an alternative 



mechanism should be invoked. Following this conclusion, 
it has been suggested that the polar caps are heated due 
to dissipation of the magnetic or/and rotational energy of 
the primary. In particular, due to the dissipation of Alfven 



waves (Wynn et al. 1997) or/and due to the impact of the 
backflowing relativistic particles accelerated in the mag- 
netosphere of the white dwarf (Ikhsanov 1998). 

A different possibility to interpret the hot pol ar cap s 
phenomenon has bee n discu ssed by Kuijpers et al. ( 1997 ), 
Meintjes & de Jager (|2000D and Choi & Yi ( |2000D . These 
authors pointed out that the mass accretion rate onto the 
white dwarf surface required to explain the observed lu- 
minosity of the pulsing component is Ma ^ lO-'^^gs"^. 
This is essentially smaller than the mass loss rate of the 
secondary (M ~ lO^^gs"^) derived from the observations 
of the Balmer continuum and the optical/UV emission 
lines (see e.g. Eracleous & Horne 1996). On this basis they 
have envisaged a situation in which the major fraction of 
the plasma inflowing into the Roche lobe of the primary 
within the orbital plane (so called low-altitude accretion 
flow) is ejected from the system due to the propeller ac- 
tion, while a small amount of material, which reaches the 
primary magnetosphere at the bases of the corotational 
cylinder (so called high-altitude accretion flow)^, is able 
to penetrate through the magnetospheric boundary and 
reach the white dwarf surface. The theoretical analysis of 
this situation (which has not been performed so far) is the 
subject of the present paper. I discuss the possible origin 
of the high-altitude accretion flow (Sect. 2), the efficiency 
of plasma penetration into the primary magnetosphere 
through the bases of the corotational cylinder (Sect. 3) 
and the trajectory of the plasma inside the white dwarf 
magnetosphere (Sect. 4). My basic conclusion is that the 
assumption about the accretion nature of the 33 s pulsing 
component in AE Aqr has no sufficient theoretical grounds 
and thus, the investigation of alternative mechanisms of 
the polar caps heating seems to be more fruitful (Sect. 5). 



2. The origin of high-altitude accretion flow 

In the general case the mass-exchange in a close binary 
system can be realized due to the stream- fed and/or 
the wind-fed mass transfer mechanisms. As shown by 



Ikhsanov ( |1997[ ) the mass capture rate by the primary 
from the stellar wind of the normal companion in AE Aqr 
is limited to Af^a ^ 3.310^^gs^^. This value is signifl- 
cantly smaller than that required to explain the observed 
UV/X-ray luminosity of the pulsing component in terms 
of accretion onto the primary surface. That is why the 
wind-fed mass transfer mechanism cannot be responsible 
for the powerful high-altitude accretion flow in AE Aqr. 

The effective cross-section of the stream flowing into 
the Roche lobe of the primary through the LI point can 



^ According to Patterson ( 1979 ) and Eracleous et al. ( 1994 ) 
the rotational axis of the white dwarf in AEAqr is almost 



^ Hereafter the distance to AE Aqr is adopted to be 100 pc perpendicular to the orbital plane 
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be limited by the cross-section of the throat at this point, 
which in the case of AEAqr is (Ikhsanov 1997) 



Q 



2ncia 



2„3 



kG{Mi +M2) 



1.85 X 10^^ 



T 



104 K 



(3) 



Here Cg and T are the sound speed and the plasma tem- 
perature at the LI point, a is the orbital separation. Mi 
and M2 are the masses of the primary and the secondary, 
respectively, and fc is a dimensionless constant depending 
on the mass ratio of the components (for discussion see 
Meyer & Meyer- Hofmeistcr 1983). This allows to evalu- 
ate the effective radius of the stream at the LI point^ as 



2.4 10^ cm 



The stream inside the Roche lobe of the primary fol- 
lows the ballistic trajectory toward the white dwarf until 
it interacts with the star magnetic field. The maximum 
possible expansion of the stream during this part of its 
trajectory is 



Ar < tffCs ~ 2 10^ Rl{ M^ 



9 ,3/2 ,,,-i/2yi/2 



0.8 



(4) 



where ts = B?/"^ / \/2GMi is the free-faU time, Rn is the 
distance from the white dwarf to the LI point expressed 
in units of 10^^ cm and T/i = T/IO"* K. Hence the altitudes 
at which the stream flows into the white dwarf magneto- 
sphere can be estimated as 



(5) 



i.e. the material inflowing through the LI point reaches 
the white dwarf magnetosphere in a narrow band around 
the rotational equator of the primary. 

The situation could be different if the secondary in 
AEAqr essentially overflows its Roche lobe. This assump- 
tion, however, contradicts the results of the investigation 
of gradual variations of optical brightness in the quiescent 
state (van Paradijs et al. 1989| ). Furthermore, in this case, 
the mass transfer rate in the system would be essentially 
higher than that currently estimated (see Introduction). 

An additional possibility to explain the origin of the 
high-altitude accretion flow is to assume that the tem- 
perature of the material surrounding the stream at the 
LI point is about of 10'' K. In this case the value of Ar 
proves to be comparable to the magnetospheric radius of 
the white dwarf and, hence, the spherically symmetrical 
accretion flow onto the primary magnetosphere can be ex- 
pected. In principle this assumption could be accepted if 
the mechanism of plasma heating at the LI point is ex- 
plained. 

3. The efficiency of plasma entry into the 
magnetosphere of the white dwarf 

As shown by Arons & Lea ( 197(^ ) the mass accretion rate 
onto the surface of a magnetized compact star is limited 



Here I would like to note that the plasma density across 
the stream at the LI point is ris oc exp{ — (r/rso)^}- 



by the rate of plasma penetration into the star magne- 
tosphere. The latter depends on the shape of the mag- 
netospheric boundary, the geometry of the accretion flow 
and the physical parameters of plasma over the boundary. 
Correspondingly, three modes of plasma penetration into 
the star magnetic field can be realized: the interchange in- 
stabilities, the diffusion, and the reconnection of the mag- 
netic field lines (for discussion see Eisner & Lamb 1984). 

Choi & Yi (2000) have distinguished two components 
of the accretion flow beyond the white dwarf magneto- 
sphere: (i) the accretion stream (the low-altitude com- 
ponent) and (ii) the almost spherically symmetrical ac- 
cretion flow (the h igh-altitude component). As shown by 
Wynn et al. (1997) the stream interacts with the magnetic 
held of the white dwarf in a local region at low altitudes 
(I K 1^ 3°) mainly at the closest approach to the pri- 
mary (rms — 10^*^ cm). The influence of the stream on the 
large scale field of the white dwarf at higher altitudes (i.e. 

o 

I K 1^ 3 ) is small and the dipole approximation for the 
field of the primary in these regions can be used. 

The interaction between the spherically symmetrical 
accretion flow and the dipole magnetic field of a com- 
pact star in the state of propeller has been investigated 



by Davies et al. (1979) and Davies & Pringle (1981). They 



have shown that the magnetosphere of the primary in this 
case is closed^ and the magnetospheric radius can be eval- 
uated as 



Tc ^ i?m = 1.8 lO^VsfA^rs^^^^W ' cm 



'1/7, 



(6) 



where /^32 = /i/lO'^^Gcm^ and Mi 5 is the mass accretion 
rate of the high-altitude component onto the magneto- 
spheric boundary expressed in units of IQi^gs^^. 

Under the condition Tc > Rear, that is just the case 
of the white dwarf in AE Aqr, the state of the primary 
is classified as supersonic propeller (the linear velocity of 
the magnetospheric boundary, — 27rrc/Ps, exceeds the 
sound speed in the accretion flow, which is limited by the 
free-fall velocity, Vs = \j2GM\lr^. In this situation the 
magnetosphere of the star is surrounded by a turbulent at- 
mosphere in which the plasma temperature is of the order 
of the free-fall temperature Tg = GMm-p/ krc (where m-p 
and k are the proton mass and the Boltzmann constant, 
respectively). 



As shown by Davies & Pringle (1981) the rate of 
plasma penetration through the magnetic field of a com- 
pact star, which is in the state of the supersonic propeller, 
is almost negligible. In the region with Tc cos k > i?cor 
(which in our case corresponds to the magnetic latitudes 

o 

I K |;^ 85 ) the centrifugal force dominates the gravi- 
tational force at the magnetospheric boundary. That is 
why the plasma entering the magnetosphere in this re- 
gion is unable to flow along the field lines to the primary 
surface but is pushed back by the fast rotating magneto- 
sphere. In the region Tc cos k ^ i?cor (which corresponds 



^ The shape of the m agnetosphere is similar to that derived 
by Arons & Lea (|l97e|) 
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to 85 ^1 K 1^ 90 ) the centrifugal force is not effec- 
tive. However, the effective penetration of plasma into the 
magnetosphere does not occur since the temperature of 
the plasma over the boundary is T(rc) ^ Tg and hence, 
the magnetospheric boundary is stable with respect to 
intercha nge in stabilities (see Arons & Lea 1976 ; Eisner 
& Lamb 1977 ). In this situation the efficiency of plasma 
penetration into the magnetic field of the primary (due 
to the magnetic field lines reconnection and/or diffusion) 
is smaller than 1% (see Ikhsanov 2001 and references 



therein). Furthermore, the area of the bases of the coro- 
tational cylinder in the considered case constitutes about 
10% of the total magnetospheric area. That is why the rate 
of plasma penetration into the white dwarf magnetosphere 
is at least three orders of magnitude smaller than the ini- 
tial mass accretion rate in the high-altitude accretion flow. 
This, however, is too small to explain the luminosity of the 
33 s pulsing component. 



4. Accretion flow inside the primary 
magnetosphere 

The last item I briefly address in this section is the plasma 
flow inside the white dwarf magnetosphere. It is well 
known that the accreting plasma penetrates the magne- 
tosphere across the filed lines forming the magnetopause, 
which is situated at the boundary and has the thickness 
<C Rm (see e.g. Arons fc Lea 1976, Ghosh & Lamb 



1979 and Eisner & Lamb 1984). Plasma penetrating the 



magnetosphere flows along the field lines within a rela- 
tively narrow channel towards the star surface. Following 
this way it can reach the star surface if all parts of the 
channel are situated inside the corotational cylinder (oth- 
erwise it will be pushed out from the magnetosphere by the 
centrifugal force). This condition is always satisfied if the 
magnetospheric radius of the star is smaller than the coro- 
tational radius, i.e. if the star is in the accrctor state. But 
if the star is in the propeller state the realization of this 
condition is not obvious. The answer depends on whether 
a magnetic line, which is situated inside the corotational 
cylinder and connects the star surface with the bases of 
the corotational cylinder, exists in the magnetosphere. 

To investigate this problem I use the dipole approxi- 
mation for the magnetic field of the primary^. A magnetic 
line of force in this case has the equation 



'A, 



(7) 



where A is the magnetic latitude and is the distance 
from the line origin to the point of its intersection with 
the plane of the magnetic equator^. I denote the angle 
between the magnetic and rotational axes of the primary 
by P and the angle between B and the radius vector by 



< r < Tc 



This approxi mation is valid in the region i?wd 
(see Arons & Lea 1976). 

® In the azimuthal direction the equation of the magnetic 
line of force is ip — const. 




Fig. 1. Schematic description of the magnetosphere of 
the white dwarf undergoing high-altitude accretion as dis- 
cussed in the text 



a. The latter can be expressed, according to Alfven 
Falthammar ( 1963| ), as 



cos A 



sma 



Vl + 3 sin^ A 



(8) 



A situation in which at least one line of force is sit- 
uated inside the corotational cylinder and connects the 
star surface with the magnetospheric boundary is shown 
in Fig. 1^. This line has a point of contact with the surface 
of the corotational cylinder in which the vector B is par- 
allel to the rotational axis of the star, ft. As it is seen from 
the figure the following conditions are satisfied in this case 



Rr 



To COS' 



A 



(9) 



(10) 



The value re for this line can be evaluated taking into 
account that it reaches the magnetospheric boundary at 
A = (7r/2-/3). Hence, 



.(A) 



cos- 



(11) 



where ^^(A) is the distance from the white dwarf to the 
base of the corotational cylinder. According to Arons & 
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Lea (1976) in the case of spherically symmetric accretion, 
^"111 (A) can be approximated as 



(cos A) 



0.2639 



for |A|<Ao, 



0.51 + 0.63(l-2iM^^^^, for | A |> Aq 



(12) 



where Aq — 80 . 

Combining Eqs. (0-|lJ) I find that plasma accretion 
onto the primary surface in the considered case could be 
realized if the following condition is satisfied 



tan A sin A 
Vl + 3sin2 A 



-(cOsA)°-2639 > 0. 



(13) 



Solving this inequality for A and putting the result to 
Eqs. (H) and ^ I find that the accretion process onto the 
surface of the white dwarf in AE Aqr could be realizedH if 

o 

/3 ^ 38 .At the same time, the value of /3 evaluated by 
Eracleous et al. (119941) from the investigation of the pulse 

' o o 

profiles of the 33 s oscillations is 75 ^ /3 ^ 77 . This 
indicates that the assumption about the direct accretion 
of plasma onto the surface of the white dwarf in AE Aqr 
contradicts the results of optical/UV observations of the 
system and thus cannot be accepted. 



5. Conclusions 

One can conclude that the assumption about the direct 
plasma accretion onto the surface of the white dwarf in 
AEAqr cannot be accepted. The theoretical grounds of 
this conclusion are the following. First, the origin of the 
high-altitude accretion flow, which moves directly to the 
bases of the corotational cylinder of the white dwarf with 
the rate > 10^^ gs~^, in the particular case of AE Aqr, is 
rather unclear. Second, even if we assume that this accre- 
tion flow exists, the efficiency of its penetration into the 
white dwarf magnetosphere is smaller than 0.1%. Finally, 
for the established value of the angle between the rota- 
tional and magnetic axes of the white dwarf in AE Aqr, 
the material penetrating the magnetosphere through the 
bases of the corotational cylinder could never reach the 
star surface flowing along the magnetic field lines. 

From the observational point of view the assumption 
about the accretion onto the surface of the white dwarf in 
AE Aqr also faces a number of serious problems. In par- 
ticular, the temperature of the plasma in the polar caps 
of the white dw arf in AEAqr (T ~ 26000 K) derived by 
Eracleous et al. (1994) is significantly smaller than the sur- 
face temperature in the magnetic pole regions of accreting 

10^ K, Eracleous et al. 



white dwarfs (typically T ^ 10** K, Eracleous et al. |l99l| ). 
Furthermore, the X-ray spectrum of AEAqr is soft and 
essentially differs from the hard X-ray spectra of all inter- 
mediate polars as well as from those of almost all accretion 



''' The corotational radius and the equatorial magnetospheric 
radius of the white dwarf are taken as _Rcor ~ 1-5 lO'' cm and 
T'c = 1.8 10^" cm, respectively (see above) 



powered close binaries (Clayton & Osborne 19951 ). On the 
other hand, due to the power law spectrum with a « —2 
and the ratio of the X-ray luminosity to the spindown lu- 
minosity Lx/Lsd ^ 10~3, the X-ray emission of AEAqr is 
rather similar to the X-rays detected from spin-powered 
pulsars (Becker & Triimper 1997). This resemblance indi- 



cates that the energy release in the magnetosphere of the 
white dwarf in AE Aqr and the magnetospheres of spin- 
powered pulsars may have a common nature and, if so, 
the assumption about the accretion onto the white dwarf 
surface proves not to be necessary. 
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